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ABSTRACT 


We present cluster counts and corresponding cosmological constraints from the Planck full mission data set. Our catalogue consists of 
439 clusters detected via their Sunyaev-Zeldovich (SZ) signal down to a signal-to-noise of six, and is more than a factor of two larger 
than the 2013 Planck cluster cosmology sample. The counts are consistent with those from 2013 and yield compatible constraints 
under the same modelling assumptions. Taking advantage of the larger catalogue, we extend our analysis to the two-dimensional 
distribution in redshift and signal-to-noise. We use mass estimates from two recent studies of gravitational lensing of background 
galaxies by Planck clusters to provide priors on the hydrostatic bias parameter, 1 - fc. In addition, we use lensing of cosmic microwave 
background (CMB) temperature fluctuations by Planck clusters as a third independent constraint on this parameter. These various 
calibrations imply constraints on the present-day amplitude of matter fluctuations in varying degrees of tension with those coming 
from Planck analysis of primary fluctuations in the CMB; for the lowest estimated values of I - b the tension is mild, only a little 
over one standard deviation, while for the largest estimated value it remains substantial. We also examine constraints on extensions to 
the base flat ACDM model by combining the cluster and CMB constraints. The combination appears to favour non-minimal neutrino 
masses, but this possibility does little to relieve the overall tension because it simultaneously lowers the implied value of the Hubble 
parameter, thereby exacerbating the discrepancy with most current astrophysical estimates. Improving the precision of cluster mass 
calibrations from the current 10%-level to 1% would significantly strengthen these combined analyses and provide a stringent test of 
the base ACDM model. 


1. Introduction 

Galaxy cluster counts are a standard cosmological tool that has 
found powerful application in recent Sunyaev-Zeldovich (SZ) 
surveys performed by the Atacama Cosmology Telescope (ACT, 
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Swetz et al. 2011; Hasselfield et al. 2013), the South Pole Tele¬ 
scope (SPT, Carlstrom et al. 2011; Benson et al. 2013; Re- 
ichardt et al. 2013; Bocquet et al. 2014), and the Planck satellite' 
(Tauber et al. 2010; Planck Collaboration I 2011). The abun- 

' Planck (http://www.esa.int/Planck) is a project of the Euro¬ 
pean Space Agency (ESA) with instruments provided by two scientific 
consortia funded by ESA member states and led by Principal Investi¬ 
gators from France and Italy, telescope reflectors provided through a 
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dance of clusters and its evolution are sensitive to the cosmic 
matter density, and the present amplitude of density fluctua¬ 
tions, characterised by erg, the rms linear overdensity in spheres 
of radius 8/1 * Mpc. The primary cosmic microwave background 
(CMB) anisotropies, on the other hand, reflect the density pertur¬ 
bation power spectrum at the time of recombination. This differ¬ 
ence is important because a comparison of the two tests the evo¬ 
lution of density perturbations from recombination until today, 
enabling us to look for possible extensions to the base ACDM 
model, such as non-minimal neutrino masses or non-zero curva¬ 
ture. 

Launched on 14 May 2009, Planck scanned the entire sky 
twice a year from 12 August 2009 to 23 October 2013 at angular 
resolutions from 33' to 5'with two instruments: the Low Fre¬ 
quency Instrument (LFI; Bersanelli et al. 2010; Mennella et al. 
2011), covering bands centred at 30, 44, and 70GHz, and the 
High Frequency Instrument (HFI; Lamarre et al. 2010; Planck 
HFI Core Team 201 1), covering bands centred at 100, 143, 217, 
353, 545, and 857 GHz. 

A initial set of cosmology results appeared in 2013 based 
on the first 15.5 months of data (Planck Collaboration I 2014), 
including cosmological constraints from the redshift distribu¬ 
tion of 189 galaxy clusters detected at signal-to-noise (SNR) 
> 7 (hereafter, our "hrst analysis" or the "2013 analysis", Planck 
Collaboration XX 2014). The present paper is part of the second 
set of cosmology results obtained from the full mission data set; 
it is based on an updated cluster sample introduced in an accom¬ 
panying paper (the PSZ2, Planck Collaboration I 2015). 

Our first analysis found fewer clusters than predicted by 
Planck’s base ACDM model, expressed as tension between the 
cluster constraints on (Qm, erg) and those from the primary CMB 
anisotropies (Planck Collaboration XVI2014). This could reflect 
the need for an extension to the base ACDM model, or indicate 
that clusters are more massive than determined by the SZ signal- 
mass scaling relation adopted in 2013. 

The cluster mass scale is the largest source of uncertainty in 
interpretation of the cluster counts. We based our first analysis 
on X-ray mass proxies that rely on the assumption of hydrostatic 
equilibrium. Simulations demonstrate that this assumption can 
be violated by bulk motions in the gas or by nonthermal sources 
of pressure (e.g., magnetic fields or cosmic rays, Nagai et al. 
2007; Piffaretti & Valdarnini 2008; Meneghetti et al. 2010). Sys- 
tematics in the X-ray analyses (e.g., instrument calibration, tem¬ 
perature structure in the gas) could also bias the mass measure¬ 
ments significantly. We quantihed our ignorance of the true mass 
scale of clusters with a mass bias parameter that was varied over 
the range [0-30]%, with a baseline value of 20% (see below for 
the definition of the mass bias), as suggested by numerical sim¬ 
ulations (see the Appendix of Planck Collaboration XX 2014). 

Gravitational lensing studies of the SZ signal-mass relation 
are particularly valuable in this context because they are inde¬ 
pendent of the dynamical state of the cluster (Marrone et al. 
2012; Planck Collaboration Int. Ill 2013), although they also, 
of course, can be affected by systematic effects (e.g., Becker & 
Kravtsov 2011). New, more precise lensing mass measurements 
for Planck clusters have appeared since our 2013 analysis (von 
der Linden et al. 2014b; Hoekstra et al. 2015). We incorporate 
these new results as prior constraints on the mass bias in the 
present analysis. Two other improvements over 2013 are use of 
a larger cluster catalogue and analysis of the counts in signal-to- 
noise as well as redshift. 


collaboration between ESA and a scientific consortium led and funded 
by Denmark, and additional contributions from NASA (USA). 



Redshift z 


Fig. 1: Mass-redshift distribution of the Planck cosmological 
samples colour-coded by their signal-to-noise, q. The baseline 
MMF3 2015 cosmological sample is shown as the small filled 
circles. Objects which were in the MMF3 2013 cosmological 
sample are marked by crosses, while those in the 2015 inter¬ 
section sample are shown as open circles. The final samples are 
defined hy q > b. The mass Myz is the Planck mass proxy (see 
text, Amaud et al. 2015). 


In addition, we apply a novel method to measure cluster 
masses through lensing of the CMB anisotropies. This method, 
presented in Melin & Bartlett (2014), enables us to use Planck 
data alone to constrain the cluster mass scale. It provides an im¬ 
portant independent mass determination, which we compare to 
the galaxy lensing results, and one that is representative in the 
sense that it averages over the entire cluster cosmology sample, 
rather than a particularly chosen subsample. 

Our conventions throughout the paper are as follows. We 
specify cluster mass, M 500 , as the total mass within a sphere 
of radius R 500 , defined as the radius within which the mean 
mass over-density of the cluster is 500 times the cosmic criti¬ 
cal density at its redshift, z: M 500 = (4;r/3)RjQQ[500pc(z)], with 
Pc{z) - 3 H^(z)I(^7tG), where H(z) is the Hubble parameter with 
present-day value Hq - hx 100 kms“^ Mpc“^. We give SZ sig¬ 
nal strength, T 500 , in terms of the Compton y-profile integrated 
within a sphere of radius R 500 , and we assume that all clusters 
follow the universal pressure profile of Arnaud et al. (2010). 
Density parameters are defined relative to the present-day crit¬ 
ical density, e.g., Qm = Pmlpdz - 0) for the total matter density. 
Pm- 

We begin in the next section with a presentation of the Planck 
2015 cluster cosmology samples. In Sect. 3 we develop our 
model for the cluster counts in both redshift and signal-to-noise, 
including a discussion of the scaling relation, scatter and the 
sample selection function. Section 4 examines the overall cluster 
mass scale in light of recent gravitational lensing measurements; 
we also present our own calibration of the cluster mass scale 
based on lensing of the CMB temperature fluctuations. Con¬ 
struction of the cluster likelihood and selection of external data 
sets is detailed in Sect. 5. We present cosmological constraints in 
Sect. 6 and then summarize and discuss our results in Sect. 7. We 
examine the potential impact of different modeling uncertainties 
in the Appendix. 
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2. The Planck cosmological samples 

We detect clusters across the six highest frequency Planck 
bands (100 - 857 GHz) using two implementations of the multi¬ 
frequency matched filter (MMF3 and MMFl, Melin et al. 2006; 
Planck Collaboration XXIX 2014) and a Bayesian extension 
(PwS, Carvalho et al. 2009) that all incorporate the known (non- 
relativistic) SZ spectral signature and a model for the spatial 
profile of the signal. The latter is taken as the so-called univer¬ 
sal pressure profile from Arnaud et al. (2010) — with the non 
standard self-similar scaling — and parameterized by an angular 
scale, 0500 . 

We empirically characterize noise (all non-SZ signals) in lo¬ 
calized sky patches (10° on a side for MMF3) using the set of 
cross-frequency power-spectra, we construct the filters with the 
resulting noise weights, and we then filter the set of six fre¬ 
quency maps over a range of cluster scales, 6500 , spanning 1 - 
35 arcmin. The filter returns an estimate of T 500 for each scale, 
based on the adopted profile template, and sources are finally as¬ 
signed the 0500 (and hence T 500 ) of the scale that maximizes their 
signal-to-noise. Details are given in Planck Collaboration XXIX 
(2014) and in an accompanying paper introducing the Planck 
full-mission SZ catalogue (PSZ2, Planck Collaboration XXVII 
2015). 

We define two cosmological samples from the general PSZ2 
catalogues, one consisting of detections by the MMF3 matched 
filter and the other of objects detected by all three methods (the 
intersection catalogue). Both are defined by a signal-to-noise 
(denoted q throughout) cut of q > b. We then apply a mask to 
remove regions of high dust emission and point sources, leaving 
65% of the sky unmasked. The general catalogues, noise maps 
and masks can be downloaded from the Planck Legacy Archive^. 

The cosmological samples can be easily constructed from 
the PSZ2 union and MMF3 catalogues. The MMF3 cosmology 
sample is the subsample of the MMF3 catalogue defined hy q> 
6 and for which the entry in the union catalogue has C0SM0=’T’. 
The intersection cosmology sample is defined from the union 
catalogue by the criteria C0SM0=’T’, PIPEDET=111, and ^ > 6. 

Fig. 1 shows the distribution of these samples in mass and 
redshift, together with the 2013 cosmology sample. The mass 
here is the Planck mass proxy, Myz, defined in Arnaud et al. 
(2015) and taken from the PSZ2 catalogue. It is calculated using 
the Planck size-flux posterior contours in conjunction with X- 
ray priors to break the size-flux degeneracy inherent to the large 
Planck beams (see, e.g. Fig. 16 of Planck Collaboration XXVII 
(2015)). The samples span masses in the range [2-10] x lO'"^ Mq 
and redshifts from z = 0 to 1^. The MMF3 (intersection) sample 
contains 439 (493) detections. Note that the intersection cata¬ 
logue has more objects than the MMF3 catalogue because of the 
different definitions of the signal-to-noise in the various cata¬ 
logues. The signal-to-noise for the intersection catalogue corre¬ 
sponds to the highest signal-to-noise of the three detection algo¬ 
rithms (MMFl, MMF3 or PwS), while for the MMF3 catalogue 
we use its corresponding signal-to-noise. As a consequence, the 
lowest value for the MMF3 signal-to-noise in the intersection 
sample is 4.8. We note that, while being above our detection 
limit, the Virgo and the Perseus clusters are not part of our sam¬ 
ples. This is because Virgo is too extended to be blindly detected 
by our algorithms and Perseus is close to a masked region. 

The 2015 MMF3 cosmology sample contains all but one of 
the 189 clusters of the 2013 MMF3 sample. The missing cluster 
is PSZl 980, which falls inside the 2015 point source mask. Six 

^ http://pla.esac.esa.int/pla/aio/planckProducts.html 
^ We fix /; = 0.7 and Ha = 1 - = 0.7 for the mass calculation. 


(14) redshifts are missing from the MMF3 (intersection) sample. 
Our analysis accounts for these by renormalizing the observed 
counts to redistribute the missing fraction uniformly across red- 
shift. The small number of clusters with missing redshifts has no 
significant impact on our results. 

We use the MMF3 cosmology sample at ^ > 6 for our base¬ 
line analysis and the intersection sample for consistency checks, 
as detailed in the Appendix. In particular, we show that the in¬ 
tersection sample yields equivalent constraints. 

3. Modelling cluster counts 

From the theoretical perspective, cluster abundance is a function 
of halo mass and redshift as specified by the mass function. Ob- 
servationally, we detect clusters in Planck through their SZ sig¬ 
nal strength or, equivalently, their signal-to-noise and measure 
their redshift with follow-up observations. The observed clus¬ 
ter counts are therefore a function of redshift, z, and signal-to- 
noise, q. While we restricted our 2013 cosmology analysis to the 
redshift distribution alone (Planck Collaboration XX 2014), the 
larger catalogue afforded by the full mission data set offers the 
possibility of an analysis in both redshift and signal-to-noise. We 
therefore develop the theory in terms of the joint distribution of 
clusters in the (z, g')-plane and then relate it to the more specific 
analysis of the redshift distribution to compare with our previous 
results. 

3.1. Counts as a function of redshift and signai-to-noise 

The distribution of clusters in redshift and and signal-to-noise 
can be written as 

dN r r dN _ 

— — I t/Gjna.sk I dAl^QQ — P\.q\qm{^5Q0^ Z, 1-, 

dzdq J J dzdMsoodLl 

( 1 ) 

with 

dN _ dN dV 
dzdMsoodQ. dVdMsoo dzd£l’ 

i.e., the dark matter halo mass function times the volume el¬ 
ement. We adopt the mass function from Tinker et al. (2008) 
throughout, apart from the Appendix where we compare to the 
Watson et al. (2013) mass function as a test of modelling robust¬ 
ness; there, we show that the Watson et al. (2013) mass function 
yields constraints similar to those from the Tinker et al. (2008) 
mass function, but shifted by about Icr towards higher and 
lower (Tg along the main degeneracy line. 

The quantity P[q\q^{Msoo,z,l,b)] is the distribution of q 
given the mean signal-to-noise value, ^tn(4^500; I, b), predicted 
by the model for a cluster of mass M 500 and redshift z lo¬ 
cated at Galactic coordinates il,b)^. This latter quantity is de¬ 
fined as the ratio of the mean SZ signal expected of a cluster, 
Y 5 ooiM 5 oo,z), as given in Eq. (7), and the detection filter noise, 
cTf [ 05 oo( 4 f 5 oo, z), I, b]: 

qm = F5oo(4^5oo,z)/O‘f[05oo(4^5oo>-Z)j>^]- (3) 

The filter noise depends on sky location {l,b) and the cluster 
angular size, 05 oo, which introduces additional dependence on 

* Note that this form assumes, as we do throughout, that the distribu¬ 
tion depends on z and M 500 only through the mean value q^ \ specifically, 
that the intrinsic scatter, (Ti„y, of Eq. (9) is constant. 
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mass and redshift. More detail on cTf can be found in Planck 
Collaboration XX (2014) (see in particular Fig. 4 therein). 

The distribution P[q\qm] incorporates noise fluctuations and 
intrinsic scatter in the actual cluster T 500 around the mean value, 
YsooiMsoo, z), predicted from the scaling relation. We discuss this 
scaling relation and our log-normal model for the intrinsic scatter 
below, and Sect. 4 examines the calibration of the overall mass 
scale for the scaling relation. 

The redshift distribution of clusters detected alq > q^^Lt is the 
integral of Eq. (1) over signal-to-noise. 


dN 

dz 


(q > ?cat) 


r“ d^ 

Xca. ^ dzdq 


/““f 


dMsooM-^soo, z, I, b) 


dN 


dzdMsoodQ. ’ 


Table 1: Summary of SZ-mass scaling law parameters (see 
Eq. 7). 


a 


b 


(4) 


Parameter 

Value 

logF. 

-0.19 ±0.02 


1.79 ±0.08 

yS*’ 

0.66 ± 0.50 

<Y 

0.127 ±0.023 


Except when specified, a is constrained by this prior in our one¬ 
dimensional likelihood over N{z), but left free in our two-dimensional 
likelihood over N{z, q). 

We fix p to its central value throughout, except when examining mod¬ 
elling uncertainties in the Appendix. 

The value is the same as in our 2013 analysis, given here in terms of 
the natural logarithm and computed from criogy = 0.075 ± 0.01. 


with 


x(M50q,z, I, b) 



dqP[q\q„,{MsQQ,Z, I, b)]. 


(5) 


Equation (4) is equivalent to the expression used in our 2013 
analysis if we write it in the form 


allows for any difference between the X-ray determined masses 
and true cluster halo mass: Mx = (1 - b)Ms(a- This is discussed 
at length in Sect. 4. 

We adopt a log-normal distribution for T 500 around its mean 
value ? 5 oo, and a delta function for ^500 centred on Ssqq: 


X ' 


J" dlnYsooJ" ^ 


dlnYsoo I d0sooP(lnY5oo,ff5Oolz,Msoo)x(Y5OO,ff5ooJ,b), 

(6) 


where xi^soo, %oo, h b) is the survey selection function at 
q > qcat in terms of true cluster parameters (Sect. 3.3), and 
/"(In T 500 , 6 * 50 ok, iWsoo) is the distribution of these parameters 
given cluster mass and redshift. We specify the relation between 
Eq. (5) and Eq. ( 6 ) in the next section. 


E(ln T 500 , 05 Oo|iW 5 OO> z) 


— _^“In^(r5oo/F5oo)/(2(T^J,Y) 

V^CrinY 

X <5[^500 - ®500]^ 


(9) 


where T 50 o(Af 5 oo,z) and 6 > 5 oo(M 5 oo,z) are given by Eqs. (7) and 
( 8 )^ The 5-function maintains the empirical definition of Rsoo 
that is used in observational determination of the the profile. 

We can now specify the relation between Eqs. (5) and ( 6 ) by 
noting that 


3.2. Observable-mass relations 


A crucial element of our modelling is the relation between clus¬ 
ter observables, T 500 and O 500 , and halo mass and redshift. Due 
to intrinsic variations in cluster properties, this relation is de¬ 
scribed by a distribution function, P(ln T 500 , ^soolA/soo, z), whose 
mean values are specified by the scaling relations Y 5 oo(Msoo,z) 
and 05 oo(I^ 5 OO,z). 

We use the same form for these scaling relations as in our 
2013 analysis: 




T>i(z)?500 

— V 

h 


m-Z7)M5oor 

10-4 Mpc2 

1 * 

0.7 


6x 1014 Mo 


(7) 


and 


^500 - 6** 


h 

- 2/3 

(1 - b)Msoo 

0.7 


3 X 1014Mo 


1/3 


Dk{z) 
500 Mpc 


, ( 8 ) 


where 6 , - 6.997 arcmin, and fiducial ranges for the parameters 
Ft, a, and p are listed in Table 1; these values are identical to 
those used in our 2013 analysis. Unless otherwise stated, we use 
Gaussian distributions with mean and standard deviation given 
by these values as prior constraints; one notable exception will 
be when we simultaneously fit for a and cosmological param¬ 
eters. In the above expressions, Df^(z) is the angular diameter 
distance and E{z) = H{z)IHq. 

These scaling relations have been established by X-ray ob¬ 
servations, as detailed in the Appendix of Planck Collaboration 
XX (2014), and rely on mass determinations, Mx, based on hy¬ 
drostatic equilibrium of the intra-cluster gas. The mass bias pa¬ 
rameter, b, assumed to be constant in both mass and redshift. 


P[q\qUM500,Z, I, b)] 


I 


d In qY[iP\^\qYa\PWn » 


( 10 ) 


where P[q\qTa\ is the distribution of observed signal-to-noise, q, 
given the model value, q^. The second distribution represents in¬ 
trinsic cluster scatter, which we write in terms of our observable- 
mass distribution, Eq. (9), as 


P[ln ^rnl^m] 


f 


i705ooP[lnT5Oo(ln^m, 


0500, 1, b), 05 Oo|iU 5 OO, z] 


V^CrinY 


( 11 ) 


Performing the integral of Eq. (5), we find 


X = 


/ 


d In q^Plln qmlqmlxiYsoo, 0500, h b). 


( 12 ) 


with the definition of our survey selection function 


X(Y500, 0500, 1, b) 


I ' 


dqP{q\qm{YsQO, 0500, U b)}. 


(13) 


We then reproduce Eq. (6) by using the first line of Eq. (11) and 
Eq. (3). 


^ In this paper, ‘In’ denotes the natural logarithm and ‘log’ the loga¬ 
rithm to base 10; the expression is written in terms of the natural loga¬ 
rithm. 
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3.3. Selection function and survey completeness 

The fundamental quantity describing the survey selection is 
P{q\qm\, introduced in Eq. (10). It gives the observed signal-to- 
noise, used to select SZ sources, as a function of model (“true”) 
cluster parameters through q'm(l 5 oo, 6 * 500 , 1 , b), and it defines the 
survey selection function xiYsoo, Osoo, l,b) via Eq. (13). We char¬ 
acterize the survey selection in two ways. The first is with an an¬ 
alytical model and the second employs a Monte Carlo extraction 
of simulated sources injected into the Planck maps. In addition, 
we perform an external validation of our selection function using 
known X-ray clusters. 

The analytical model assumes pure Gaussian noise, in which 
case we simply have P[q\qm\ - V^. The survey se¬ 

lection function is then given by the Error Eunction (ERE), 


^( 7500 ,^ 500 ,^,^) 


1 - erf 


gcat - ^111(^500, ^ 500 , 1 , b ) 

V2 


(14) 


This model can be applied to a catalogue with well-defined noise 
properties, i.e., cTf, such as our MME3 catalogue, but not to the 
intersection catalogue based on the simultaneous detection with 
three different methods. This is our motivation for choosing the 
MME3 catalogue as our baseline. 

In the Monte Carlo approach, we inject simulated clusters 
directly into the Planck maps and (re)extract them with the com¬ 
plete detection pipeline. Details are given in the accompanying 
2015 SZ catalogue paper, Planck Collaboration XXVII (2015). 
This method provides a more comprehensive description of the 
survey selection by accounting for a variety of effects beyond 
noise. In particular, we vary the shape of the SZ profile at fixed 
T 500 and 0500 to quantify its effect on catalogue completeness. 

We also perform an external check of the survey complete¬ 
ness using known X-ray clusters from the MCXC compilation 
(Piffaretti et al. 2011) and also SPT clusters from Bleem et al. 
(2014). Details are given in the 2015 SZ catalogue paper, Planck 
Collaboration XXVII (2015). Eor the MCXC compilation, we 
rely on the expectation that at redshifts z < 0.2 any Planck- 
detected cluster should be found in one of the ROS AT catalogues 
(Chamballu et al. 2012), because at low redshift ROSAT probes 
to lower masses than Planck^. The MCXC catalogue provides 
a truth table, replacing the input cluster list of the simulations, 
and we compute completeness as the ratio of objects in the cos¬ 
mology catalogue to the total number of clusters. As discussed 
in Planck Collaboration XXVII (2015), the results are consis¬ 
tent with Gaussian noise and bound the possible effect of profile 
variations. We arrive at the same conclusion when applying the 
technique to the SPT catalogue. 

Planck Collaboration XXVII (2015) discusses completeness 
checks in greater detail. One possible source of bias is the 
presence of correlated IR emission from cluster member galax¬ 
ies. Planck Collaboration XXIII (2015) suggests that IR point 
sources may contribute significantly to the cluster SED at the 
Planck frequencies, especially at higher redshift. The potential 
impact of this effect warrants further study in future work. 

This gives us different estimations of the selection function 
for MMP3 and the intersection catalogues. We test the sensitiv¬ 
ity of our cosmological constraints to the selection function in 
the Appendix by comparing results obtained with the different 
methods and catalogues. We find that our results are insensitive 
to the choice of completeness model (Pig A. 1), and we therefore 

® In fact, this expectation is violated to a small degree. As discussed 
in Planck Collaboration XXVII (2015), there appears to be a small pop¬ 
ulation of X-ray under-luminous clusters. 


adopt the analytical ERE completeness function for simplicity 
throughout the paper. 


4. The cluster mass scale 

The characteristic mass scale of our cluster sample is the critical 
element in our analysis of the counts. It is controlled by the mass 
bias factor, I - b, accounting for any difference between the X- 
ray mass proxies used to establish the scaling relations and the 
true (halo) mass: Mx - (1 - 6 >)M 5 oo. Such a difference could arise 
from cluster physics, such as a violation of hydrostatic equilib¬ 
rium or temperature structure in the gas, from observational ef¬ 
fects, e.g., instrumental calibration, or from selection effects bi¬ 
asing the X-ray samples relative to SZ- or mass-selected samples 
(Angulo et al. 2012). 

In our 2013 analysis, we adopted a flat prior on the mass bias 
over the range I - b — [0.7,1.0], with a reference model defined 
by 1 -b = 0.8. This was motivated by a comparison of the Y-Mx 
relation with published Y - M relations derived from numerical 
simulations, as detailed in the Appendix of Planck Collaboration 
XX (2014); this estimate was consistent with most predictions 
for any violation of hydrostatic equilibrium (although not all) 
as well as observational constraints from the available lensing 
observations. Effects other than cluster physics can contribute to 
the mass bias, as discussed in the paper, and as emphasized by 
the survey of cluster multi-band scaling relations by Rozo et al. 
(2014a,b,c). 

The mass bias was the largest uncertainty in our 2013 anal¬ 
ysis, and it severely hampered understanding of the tension 
found between constraints from the primary CMB and the cluster 
counts. Here, we incorporate new lensing mass determinations of 
Planck clusters to constrain the mass bias. We also apply a novel 
method to measure object mass based on lensing of CMB tem¬ 
perature anisotropies behind clusters (Melin & Bartlett 2014). 
These constraints are used as prior information in our analysis 
of the counts. As we will see, however, uncertainty in the mass 
bias remains our largest source of uncertainty, mainly because 
these various determinations continue to differ from 10 to 30%. 

In general, the mass bias could depend on cluster mass and 
redshift, although we will model it by a constant in the follow¬ 
ing. Our motivation is one of practicality: the limited size and 
precision of current lensing samples makes it difficult to con¬ 
strain any more than a constant, i.e., the overall mass scale of 
our catalogue. Large lensing surveys like Euclid, WEIRST, and 
the Large Synoptic Survey Telescope, and CMB lensing will im¬ 
prove this situation in coming years. 

4.1. Constraints from gravitational shear 

Several cluster samples with high quality gravitational shear 
mass measurements have appeared since 2013. Among these, 
the Weighing the Giants (WtG, von der Linden et al. 2014a), 
CLASH (Postman et al. 2012; Merten et al. 2014; Umetsu et al. 
2014), and the Canadian Cluster Comparison Project (CCCP, 
Hoekstra et al. 2015) programmes offer constraints on our mass 
bias factor, 1-b, through direct comparison of the lensing masses 
to the Planck mass proxy, Myz- 

The analysis by the WtG programme of 22 clusters from the 
2013 Planck cosmology sample yields I - b - 0.688 + 0.072. 
Their result lies at the very extreme of the range explored in 
Planck Collaboration XX (2014) and would substantially re¬ 
duce the tension found between primary CMB and galaxy clus¬ 
ter constraints. Hoekstra et al. (2015) report a smaller bias of 
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Table 2: Summary of mass scale priors 



Prior name 

Quantity 

Value & Gaussian errors 

Weighing the Giants (WtG) 
Canadian Cluster Comparison 

l-b 

0.688 ± 0.072 

Project (CCCP) 

\-b 

0.780 ± 0.092 

CMB lensing (LENS) 

1 /(1-z.) 

0.99 ±0.19 

Baseline 2013 

\-b 

0.8 [-0.1,+0.2] 


Notes. CMB lensing directly measures 1/(1 - b), which we implement 
in our analysis; purely for reference, that constraint translates approxi¬ 
mately to 1 - Zi = l.Ol/jj ij. The last line shows the 2013 baseline — a 
reference model defined by 1 - Z> = 0.8 with a flat prior in the [0.7, 1] 
range. 


Fig. 2; The cluster mass scale determined by CMB lensing. We 
show the ratio of cluster lensing mass. Miens, to the SZ mass 
proxy, Myz, as a function of the mass proxy for clusters in the 
MMF3 2015 cosmology sample. The cluster mass is measured 
through lensing of CMB temperature anisotropies in the Planck 
data (Melin & Bartlett 2014). Individual mass measurements 
have low signal-to-noise, but we determine a mean ratio for the 
sample of Miens/Myz - Ijil - b) - 0.99 + 0.19. For clarity, only 
a fraction of the error bars are plotted (see text). 


\-b - 0.78+0.07 (stat) +0.06 (sys) for a set of 20 common clus¬ 
ters, which is in good agreement with the fiducial value adopted 
in our 2013 analysis. In our analysis we add the statistical and 
systematic uncertainties in quadrature (see Table 2). 

4.2. Constraints from CMB lensing 

Measuring cluster mass through CMB lensing has been dis¬ 
cussed in the literature for some time since the study performed 
by Zaldarriaga & Seljak (1999) (see also Lewis & Challinor 
2006). We apply a new technique for measuring cluster masses 
through lensing of CMB temperature anisotropies (Melin & 
Bartlett 2014), allowing us to calibrate the scaling relations us¬ 
ing only Planck data. This is a valuable alternative to the galaxy 
lensing observations because it is independent and affected by 
different possible systematics. Additionally, we can apply it to 
the entire cluster sample to obtain a mass calibration represen¬ 
tative of an SZ flux selected sample. Similar approaches using 
CMB lensing to measure halo masses were recently applied by 
SPT (Baxter et al. 2014) and ACT (Madhavacheril et al. 2014). 

Our method first extracts a clean CMB temperature map with 
a constrained internal linear combination (ILC) of the Planck 
frequency channels in the region around each cluster; the ILC is 
constrained to nullify the SZ signal from the clusters themselves 
and provide a clean CMB map of 5 arcmin resolution. Using a 
quadratic estimator on the CMB map, we reconstruct the lensing 
potential in the field and then filter it to obtain an estimate of the 
cluster mass. The filter is an NFW profile (Navarro et al. 1997) 
with scale radius set by the Planck mass proxy for each clus¬ 
ter, and designed to return an estimate of the ratio Miens/Myz, 
where Myz is the Planck SZ mass proxy. These individual mea¬ 
surements are corrected for any mean-field bias by subtracting 
identical filter measurements on blank fields; this accounts for 
effects of apodization over the cluster fields and correlated noise. 
The technique has been tested on realistic simulations of Planck 
frequency maps. More detail can be found in Melin & Bartlett 
(2014). 


Figure 2 shows Miens/Myz as a function of Myz for all clus¬ 
ters in the MMF3 cosmology sample. Each point is an individual 
cluster^. For clarity, only some of the error bars on the ratio are 
shown; the error bars vary from 1.8 at the high mass end to 8.5 
at the low mass end with a median of 4.2. There is no indication 
of a correlation between the ratio and Myz, and we therefore fit 
for a constant ratio of Miens/Myz by taking the weighted mean 
(using the individual measurement uncertainties as provided by 
the filter) over the full data set. If the ratio differs from unity, we 
apply a correction to account for the fact that our filter aperture 
was not perfectly matched to the clusters. The correction is cal¬ 
culated assuming an NFW profile and is the order of a percent. 

The final result is 1/(1 - h) = 0.99 + 0.19, traced by the 
blue band in the figure. Note that the method constrains 1/(1 - 
b) rather than 1 - h as in the case of the shear measurements. 
The calculated uncertainty on the weighted mean is consistent 
with a bootstrap analysis where we create new catalogues of the 
same size as the original by sampling objects from the original 
catalogue with replacement; the uncertainty from the bootstrap 
is then taken as the standard deviation of the bootstrap means. 

4.3. Summary 

The three mass bias priors are summarized in Table 2, and we 
will compare cosmological constraints obtained from each. We 
will assume Gaussian distributions for 1 - h (gravitational shear) 
or 1/(1 - b) (CMB lensing) with standard deviations given by 
the error column. We favour these three lensing results because 
of their direct comparison to the Planck mass proxy. 

5. Analysis methodology 

5.1. Likelihood 

Our 2013 analysis employed a likelihood built on the cluster red- 
shift distribution, dNjdz. With the larger 2015 catalogue, our 
baseline likelihood is now constmcted on counts in the {z,q)- 
plane. We divide the catalogue into bins of size Az = 0.1 (10 
bins) and A log q - 0.25 (5 bins), each with an observed number 
N{Zi, qj) - Nij of clusters. Modelling the observed counts, Njj, 
as independent Poisson random variables, our log-likelihood is 

N,N, 

InL = 2 K'lnMv - Nij - ln[/V,/]], (15) 

ij 


’’ The values can be negative due to noise fluctuations and the low 
signal-to-noise of the individual measurements. 
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where and Nq are the total number of redshift and signal-to- 
noise bins, respectively. The mean number of objects in each bin 
is predicted by theory according to Eq. (1): 

dN 

Nij = ■^^(zi,qj)AzAq, (16) 

which depends on the cosmological (and cluster modelling) 
parameters. In practice, we use a Monte Carlo Markov chain 
(MCMC) to map the likelihood surface around the maximum 
and establish confidence limits. 

Eq. (15) assumes the bins are uncorrelated, while a more 
complete description would include correlations due to large- 
scale clustering. In practice, our cluster sample contains mostly 
high mass systems for which the impact of these effects is weak 
(e.g., Hu & Kravtsov 2003, in particular their Fig. 4 for the im¬ 
pact on constraints in the plane). 

5.2. External data sets 

Cluster counts cannot constrain all pertinent cosmological pa¬ 
rameters. They are most sensitive to Qm and erg, and when 
analysing the counts alone we must apply additional observa¬ 
tional constraints as priors on other parameters. For this pur¬ 
pose, we adopt Big Bang nucleosynthesis (BBN) constraints 
from Steigman (2008), Q.i,h^ - 0.022 + 0.002, and constraints 
from baryon acoustic oscillations (BAO). The latter combine the 
6dF Galaxy Survey (Beutler et al. 201 1), the SDSS Main Galaxy 
Sample (Padmanabhan et al. 2012; Anderson et al. 2012) and 
the BOSS DRll (Anderson et al. 2014). We refer the reader to 
Sect. 5.2 in Planck Collaboration XIII (2015) for details of the 
combination. We also include a prior on its from Planck Collab¬ 
oration XVI (2014), ris = 0.9624 ± 0.014. When explicitly spec¬ 
ified in the text, we add the supernovas constraint from SNLS- 
II and SNLS3: the Joint Light-curve Analysis constraint (JLA, 
Betoule et al. 2014). The BAO are particularly sensitive to Hq, 
while the supernovas allow precise constraints on the dark energy 
equation-of-state parameter, w. 


6. Cosmological constraints 

We begin by verifying consistency with the results of Planck 
Collaboration XX (2014) (Sect. 6.1) based on the one¬ 
dimensional likelihood over the redshift distribution, dNjdz 
(Eq. 4). We then examine the effect of changing to the full two- 
dimensional likelihood, dNjdzdq (Eq. 1) in Sect. 6.2. For this 
purpose we compare constraints on the total matter density, 
and the linear-theory amplitude of the density perturbations to¬ 
day, cTg, using the cluster counts in combination with external 
data and fixing the mass bias. 

The two-dimensional likelihood dN I dzdq is then adopted as 
the baseline in the rest of the paper. We extract constraints on 
Qm and erg from the cluster counts in combination with exter¬ 
nal data, imposing the different cluster mass scale calibrations 
as prior distributions on the mass bias. After comparing these 
new constraints to those from the CMB anisotropies in the base 
ACDM model (Sect. 6.3), we move to joint analysis of the clus¬ 
ter counts and CMB anisotropies to study, for example, exten¬ 
sions to the base ACDM model, such as non-minimal neutrino 
mass (Sect. 6.4). In these studies we vary all six parameters of 
the (flat) base ACDM model, except when considering model ex¬ 
tensions for which we include the relevant parameters. Figures 
display contours delineating marginalized constraints. 


6.1. Constraints on and erg; comparison to 2013 

Figure 3 presents constraints from the MMF3 cluster counts 
combined with the BAO and BBN priors of Sect. 5.2; we re¬ 
fer to this data combination as “SZ-i-BAO-i-BBN”. To compare 
to results from our 2013 analysis (the grey, filled ellipses), we 
use a one-dimensional likelihood based on Eq. (4) over the red¬ 
shift distribution and have adopted the reference scaling relation 
of 2013, i.e., Eqs. (7) and (8) with I - b = 0.8. For the present 
comparison, we use the updated BAO constraints discussed in 
Sect. 5.2; these are stronger than the BAO constraints used in 
the 2013 analysis, and the grey contours shown here are conse¬ 
quently smaller than in Planck Collaboration XX (2014). 

Limiting the 2015 catalogue to q > 8.5 produces a sam¬ 
ple with 190 clusters, similar to the 2013 cosmology catalogue 
(189 objects). The two sets of constraints demonstrate good con¬ 
sistency, and they remain consistent while becoming tighter as 
we decrease the signal-to-noise threshold of the 2015 catalogue. 
Under similar assumptions, our 2015 analysis thus confirms the 
2013 results reported in Planck Collaboration XX (2014). 

The area of the ellipse from ^ = 8.5 to g' = 6 decreases by 
a factor of 1.3. This is substantially less than the factor of 2.3 
expected from the ratio of the number of objects in the two sam¬ 
ples. This difference may be related to the decreasing goodness- 
of-fit of the best model as the signal-to-noise decreases. 

Figure 4 overlays the observed cluster redshift distribution 
on the predictions from the best-fit model in each case. We see 
that the models do not match the counts in the second and third 
redshift bins (counting from z = 0), and that the discrepancy, 
already marginally present at the high signal-to-noise cut corre¬ 
sponding to the 2013 catalogue, becomes more pronounced to¬ 
wards the lower signal-to-noise thresholds. This dependence on 
signal-to-noise may suggest that the data prefer a different slope, 
a, of the scaling relation than allowed by the prior of Table 1. 
We explore the effect of relaxing the X-ray prior on a in the next 
section. 

6.2. Constraints on and erg: two-dimensional analysis 

In Fig. 5 we compare constraints from the one- and two- 
dimensional likelihood with a either free or with the prior 
of Table 1. For this comparison, we continue with the 
“SZh-BAO-hBBN” data set, but adopt the CCCP prior for the 
mass bias and only consider the full 2015 MMF3 catalogue at 
q > 6. 

The grey and black contours and lines show results from the 
one-dimensional likelihood fit to the redshift distribution using, 
respectively, the X-ray prior on a and leaving a free. The redshift 
counts do indeed favour a steeper slope, and we find a posterior 
of O' = 2.28 + 0.17 in the latter case. 

We define a generalized measure of goodness-of-fit as 
~ ’ determining the probability to exceed 

{PTE) the observed value using Monte Carlo simulations of 
Poisson statistics for each bin with the best-fit model mean A,. 
The observed value of the fit drops from 17 {PTE - 0.07) with 
the X-ray prior, to 15 {PTE - 0.11) when leaving a free. When 
leaving a free, Gm increases and erg decreases, following their 
correlation with a shown by the contours, and their uncertainty 
increases due to the added parameter. 

The two-dimensional likelihood over dNjdzdq better con¬ 
strains the slope when a is free, as shown by the violet curves 
and contours. In this case, the preferred value drops back towards 
the X-ray prior: a = 1.89 + 0.12, just over 1 cr from the central 
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Fig. 3; Contours at 2cr for different signal-to-noise thresholds, 
q - 8.5, 7, and 6, applied to the 2015 MMF3 cosmology sam¬ 
ple for the SZ-hBAO-hBBN data set. The contours are compati¬ 
ble with the 2013 constraints (Planck Collaboration XX 2014), 
shown as the filled, light grey ellipses at 1 and 2cr (for the BAO 
and BBN priors of Sect 5.2; see text). The 2015 catalogue thresh- 
olded 2 it q > 8.5 has a similar number of clusters (190) as the 
2013 catalogue (189). This comparison is made using the an¬ 
alytical etTor-function model for completeness and adopts the 
reference observable-mass scaling relation of the 2013 analysis 
- b - 0.8, see text]. The redshift distributions of the best-fit 
models are shown in Fig. 4. For this figure and Fig. 4, we use the 
one-dimensional likelihood over the redshift distribution, dNjdz 
(Eq. 4). 


X-ray value. Re-imposing the X-ray prior on a with the two- 
dimensional likelihood (blue curves) does little to change the 
parameter constraints. Although the one-dimensional likelihood 
prefers a steeper slope than the X-ray prior, the two-dimensional 
analysis does not, and the cosmological constraints remain ro¬ 
bust to varying a. 

We define a generalized statistic as described above, now 
over the two-dimensional bins in the (z, q')-plane. This general¬ 
ized x^ for the fit with the X-ray prior is 43 (PTE - 0.28), com¬ 
pared tox^ - 45 (PTE - 0.23) when a is a free parameter. 

Fig. 6 displays the redshift distribution of the best-fit models 
in all four cases. Despite their apparent difficulty in matching the 
second and third redshift bins, the PTE values suggest that these 
fits are moderately good to acceptable. Note that, as mentioned 
briefly in Sect. 5.1, clustering effects will increase the scatter 
in each bin slightly over the Poisson value we have assumed, 
causing our quoted PTE values to be somewhat smaller than the 
true ones. 

6.3. Constraints on and cr^: comparison to primary CMB 

Our 2013 analysis brought to light tension between constraints 
on Qm and erg from the cluster counts and those from the primary 
CMB in the base ACDM model. In that analysis, we adopted a 
flat prior on the mass bias over the range 1 - b = [0.7,1.0], with 
a reference model defined hy \ - b - 0.8 (see discussion in the 
Appendix of Planck Collaboration XX 2014). Given the good 
consistency between the 2013 and 2015 cluster results (Pig. 3), 
we expect the tension to remain under the same assumptions con¬ 
cerning the mass bias. 



0.0 0.2 0.4 0.6 0.8 1.0 


z 

Pig. 4: Comparison of observed counts (points with etTor bars) 
with predictions of the best-fit models (solid lines) from the 
one-dimensional likelihood for three different thresholds applied 
to the 2015 MMP3 cosmology sample. The mismatch between 
observed and predicted counts in the second and third lowest 
redshift bins, already noticed in the 2013 analysis, increases at 
lower thresholds, q. The best-fit models are defined by the con¬ 
straints shown in Pig. 3. Por this figure and Pig. 3, we use our 
one-dimensional likelihood over the redshift distribution, dNjdz 
(Eq. 4). 



Pig. 5; Comparison of constraints from the one-dimensional 
(dNjdz) and two-dimensional (dNjdzdq) likelihoods on cosmo¬ 
logical parameters and the scaling relation mass exponent, a. Por 
this comparison, we adopt the CCCP prior on the mass bias and 
the SZ-i-BAO-i-BBN data set. The corresponding best-fit model 
redshift distributions are shown in Pig. 6. 


Pigure 7 compares our 2015 cluster constraints (MMP3 
SZ-hBAO-hBBN) to those for the base ACDM model from the 
Planck CMB anisotropies. The cluster constraints, given the 
three different priors on the mass bias, are shown by the filled 
contours at 1 and 2cr, while the dashed black contours give the 
Planck TT, TP, EE-i-lowP constraints (hereafter Planck primary 
CMB, Planck Collaboration XIII 2015); the grey shaded regions 
add BAO to the CMB. The central value of the WtG mass prior 
lies at the extreme end of the range used in 2013 (i.e., l-b=0.7); 
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Fig. 6: Redshift distribution of best-fit models from the four 
analysis cases shown in Fig. 5. The observed counts in the 
MMF3 catalogue (q > 6) are plotted as the red points with error 
bars, and as in Fig. 5 we adopt the CCCP mass prior with the 
SZ-i-BAOh-BBN data set. 
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Fig. 7: Comparison of constraints from the CMB to those from 
the cluster counts in the (Qm, (T8)-plane. The green, blue and 
violet contours give the cluster constraints (two-dimensional 
likelihood) at 1 and 2cr for the WtG, CCCP, and CMB lens- 
ing mass calibrations, respectively, as listed in Table 2. These 
constraints are obtained from the MMF3 catalogue with the 
SZ-i-BAOh-BBN data set and a free. Constraints from the Planck 
TT, TE, EE-hIowP CMB likelihood (hereafter, Planck primary 
CMB) are shown as the dashed contours enclosing 1 and 2cr con- 
hdence regions (Planck Collaboration XIII2015), while the grey 
shaded region also include BAO. The red contours give results 
from a joint analysis of the cluster counts, primary CMB and 
the Planck lensing power spectrum (Planck Collaboration XV 
2015), leaving the mass bias parameter free and a constrained 
by the X-ray prior. 



with its uncertainty range extending even lower, the tension with 
primary CMB is greatly reduced, as pointed out by von der Lin¬ 
den et al. (2014b). With similar uncertainty but a central value 
shifted to 1 - h = 0.78, the CCCP mass prior results in greater 
tension with the primary CMB. The lensing mass prior, finally, 
implies little bias and hence much greater tension. 


Pig. 8: Comparison of cluster and primary CMB constraints in 
the base ACDM model expressed in terms of the mass bias, 
I - b. The solid black curve shows the distribution of values re¬ 
quired to reconcile the counts and primary CMB in ACDM; it 
is found as the posterior on the \ — b from a joint analysis of 
the Planck cluster counts and primary CMB when leaving the 
mass bias free. The coloured dashed curves show the three prior 
distributions on the mass bias listed in Tab. 2. 


6.4. Joint Planck 2014 primary CMB and cluster constraints 

We now turn to a joint analysis of the cluster counts and primary 
CMB. We begin by hnding the mass bias required to remove ten¬ 
sion with the primary CMB, and then consider one-parameter 
extensions to the base ACDM model, varying the curvature, the 
Thomson optical depth to reionization, the dark energy equation- 
of-state, and the neutrino mass scale. Unless otherwise stated, 
"CMB" in the following means Planck TT, TE, EE-i-lowP as de- 
hned in Planck Collaboration XIII (2015). All intervals are 68% 
conhdence and all upper/lower limits are 95%. 

6.4.1. Mass bias required by CMB 

In Pig. 8 we compare the three prior distributions to the mass 
bias required by the primary CMB. The latter is obtained as the 
posterior on (1 — b) from a joint analysis of the MMP3 cluster 
counts and the CMB with the mass bias as a free parameter. The 
best-fit value in this case is (1 - h) = 0.58 + 0.04, more than Icr 
below the central WtG value. Perfect agreement with the primary 
CMB would imply that clusters are even more massive than the 
WtG calibration. This figure most clearly quantihes the tension 
between the Planck cluster counts and primary CMB. 

6.4.2. Curvature 

By itself the CMB only poorly determines the spatial curvature 
(Sect. 6.2.4 of Planck Collaboration XIII 2015), but by including 
another astrophysical observation, such as cluster counts, it can 
be tightly constrained. Our joint cluster and CMB analysis, with¬ 
out external data, yields Ok = -0.012 + 0.008, consistent with 
the constraint from Planck CMB and BAO Ok = 0.000 + 0.002. 

6.4.3. Reionization optical depth 

Primary CMB temperature anisotropies also provide a precise 
measurement of the parameter combination where t is 
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Fig. 9: Constraints on the reionization optical depth, t. The 
dashed black curve is the constraint from Planck CMB (i.e. 
TT, TE, EE+lowP), while the three coloured lines are the pos¬ 
terior distribution on r from a joint analysis of the cluster counts 
and Planck TT only for the three different mass bias parameters. 


the optical depth from Thomson scatter after reionization and Ag 
is the power spectrum normalization on large scales (Planck Col¬ 
laboration XIII 2015). Low-^ polarization anisotropies break the 
degeneracy by constraining t, but this measurement is delicate 
given the low signal amplitude and difficult systematic effects; it 
is important, however, in the determination of erg. It is therefore 
interesting to compare the Planck primary CMB constraints on 
T to those from a joint analysis of the cluster counts and primary 
CMB without the low-^ polarization data (lowP). Battye et al. 
(2014), for instance, pointed out that a lower value for r than 
suggested by WMAP could reduce the level of tension between 
CMB and large scale structure. 

The comparison is shown in Pig. 9. We see that the Planck 
TT -I- SZ constraints are in good agreement with the value from 
Planck CMB (i.e., TT, TE, EE-i-lowP), with the preferred value 
for WtG slightly higher and CMB lensing pushing towards a 
lower value. The ordering CMB lensing/CCCP/WtG from lower 
to higher t posterior values matches the decreasing level of ten¬ 
sion with the primary CMB on erg. These values remain, how¬ 
ever, larger than what is required to fully remove the tension 
in each case. The posterior distributions for the mass bias are 
l-b^ 0.60 + 0.042, l-b = 0.63 + 0.047, l-b^ 0.66 + 0.045, 
respectively, for WtG, CCCP and CMB lensing, all significantly 
shifted from the corresponding priors of Table 2. Allowing r to 
adjust offers only minor improvement in the tension reflected 
by Pig. 8. Interestingly, the Planck TT posterior shown in Pig. 8 
of Planck Collaboration XIII (2015) peaks at significantly higher 
values, while our Planck TT -hSZ constraints are consistent with 
the constraint from Planck TT + lensing, an independent con¬ 
straint on T without lowP. 

6.4.4. Dark energy 

In Pig. 10 we examine constraints on a constant dark energy 
equation-of-state parameter, w. Analysis of the primary CMB 
alone results in highly degenerate grey contours. The degener¬ 
acy is broken by adding constraints such as BAG (blue contours) 
or supernovae distances (light blue contours), both picking val¬ 
ues around w = -1. The SZ counts (two-dimensional likelihood 
with CCCP prior) only marginally break the degeneracy when 
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Pig. 10; Constraints on a constant dark energy equation-of-state 
parameter, w. Analysis of the primary CMB alone yields the grey 
contours that are highly degenerate. Adding either BAG or su¬ 
pernovae to the CMB breaks the degeneracy, giving constraints 
around w = -1. Adding SZ cluster counts from the MMP3 2015 
catalogue instead to the CMB results in the rose-coloured con¬ 
tours. The green contours are constraints from joint analysis of 
the SZ counts and BAG; although much less constraining they 
agree with the CMB-i-JLA combinations and are completely in¬ 
dependent. 


combined with the CMB, but when combined with BAG they 
do yield interesting constraints (green contours) that are con¬ 
sistent with the independent constraints from the primary CMB 
combined with supernovae. We obtain Q,„ = 0.313 + 0.025 and 
w = -1.00 + 0.18 for SZ-hBAG, and Q,„ = 0.306 + 0.013 and 
w = -1.10 + 0.06 for CMB+BAG. 

6.4.5. 2 lUy 

An important, well-motivated extension to the base ACDM 
model that clusters can help constrain is a non-minimal sum 
of neutrino masses, 2 niy > 0.06 eV. Given the primary CMB 
anisotropies, the amplitude of the density perturbations today, 
characterized by the equivalent linear theory extrapolation, cg, 
is model dependent; it is a derived parameter depending, for ex¬ 
ample, on the composition of the matter content of the universe. 
Cluster abundance, on the other hand, provides a direct measure¬ 
ment of CTg at low redshifts, and comparison to the value derived 
from the CMB tests the adopted cosmological model. 

By free-streaming, neutrinos damp the growth of matter per¬ 
turbations. Gur discussion thus far has assumed the minimum 
mass for the three known neutrino species. Increasing their mass, 
2 > 0.06 eV, lowers erg because the neutrinos have larger 

gravitational influence on the total matter perturbations. This 
goes in the direction of reconciling tension — the strength of 
which depends on the mass bias — between the cluster and pri¬ 
mary CMB constraints. Cluster abundance, or any measure of erg 
at low redshift, is therefore an important cosmological constraint 
to be combined with those from the primary CMB. 

Figure 11 presents a joint analysis of the cluster counts for 
the CCCP mass bias prior with primary CMB, the Planck lens¬ 
ing power spectrum and BAG. The results without BAG (green 
and red shaded contours) allow relatively large neutrino masses, 
up to 2 ~ 0.5 eV; and when adding the lensing power spec¬ 

trum, a small, broad peak appears in the posterior distribution 
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Fig. 11: Parameter constraints on the ACDM+non-minimal neutrino mass model. For this study, we adopt the CCCP prior on the 
mass bias (see Tab. 2) and leave the scaling exponent, a, free. The green and red shaded regions show, respectively, the 1 and 2cr 
confidence regions for joint analyses of the cluster counts using the primary CMB, and the primary CMB plus the lensing power 
spectrum. The solid and dashed black contours add to these two cases constraints from BAO. 


just above ^ = 0.2 eV. We also notice some interesting cor¬ 

relations; the amplitude, erg, anti-correlates with neutrino mass, 
as does the Hubble parameter, and larger values of a correspond 
to larger neutrino mass, lower Nq and lower erg. 

As discussed in detail in Planck Collaboration XIII (2015), 
the anti-correlation with the Hubble parameter maintains the ob¬ 
served acoustic peak scale in the primary CMB. Increasing neu¬ 
trino mass to simultaneously accommodate the cluster and pri¬ 
mary CMB constraints by lowering erg, while allowed in this 
joint analysis, would therefore necessarily increase tension with 
some direct measurements of No (see discussion in Planck Col¬ 
laboration XIII (2015)). Including the BAO data greatly restricts 
this possibility, as shown by the solid and dashed black curves. 

The solid and dashed, red and black curves in Fig. 12 re¬ 
produce the marginalized posterior distributions on 2 from 
Fig. 11 . The solid blue curve is the result of a similar analysis 
where in addition the artificial parameter Al is allowed to vary. 
This parameter characterizes the amount of lensing in the tem¬ 
perature power spectrum relative to the best fit model (Planck 
Collaboration XIII 2015). Planck TT + lowP alone constraints 

Al = 1.22 + 0.10 

in mild tension with the value predicted for the ACDM model, 
Al = 1. In the base ACDM model, this parameter is fixed to 
unity, but it is important to note it is degenerate with ^ Left 
free, it allows less lensing power, which is also in line with the 
direct measurement of the lensing power spectrum (labelled as 


Lensing PS) from the four-point function (see Planck Collab¬ 
oration XIII 2015). In that light, we see that adding Al as a 
free parameter accentuates the peak in the CMB+SZ+Lensing 
PS posterior. The small internal tension between CMB+SZ and 
CMB+SZ+Al posteriors may point towards a need for an exten¬ 
sion of the minimal six-parameter ACDM. 

These posteriors lead to the following constraints: Y^niy < 
0.53 eV (95%) for CMB+SZ+Lensing PS and < 0.22 eV 
(95%) for CMB+SZ+BAO. 

We may compare these with the constraints from the primary 
CMB presented in Planck Collaboration XIII (2015). The Planck 
primary CMB by itself places an upper limit of Z '^v < 0-49 eV 
(95%), and the addition of BAO tightens this to Z '«v < 0.17 eV 
(95%). Addition of the Planck lensing power spectrum to the 
primary CMB weakens the constraint to 2 my < 0.59 eV (95%), 
as we would expect given the results and discussion above. The 
final constraint adopted by Planck Collaboration XIII (2015), for 
its robustness to possible remaining low level systematics in the 
polarization data, is 2 my < 0.23 eV (95%), not too different 
from the peak suggested in CMB+SZ+lensing PS posterior. 

Adding neutrino mass should lower erg, letting it move to¬ 
wards values favoured by the cluster counts. We might expect 
that the CMB+SZ combination would therefore find clear evi¬ 
dence for non-minimal neutrino mass. In spite of this, the green 
curve only places an upper limit on Z ^v- We may understand 
this by looking at the posterior on the mass bias I - b. The 
allowed values are well separated from the prior distribution. 
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Fig. 12: Constraints on 2 niy from a joint analysis of the clus¬ 
ter counts and primary CMB. The solid and dashed, red and 
black lines reproduce the marginalized posterior distributions 
from Fig. 11 . The solid blue line is the posterior from a similar 
analysis, but marginalized over the additional parameter Al (see 
text). If applied to the present Planck cluster cosmology sample, 
a future mass calibration of I - b = 0.80 + 0.01 would result in 
the bold, dotted black posterior curve. 


Fig. 13: Prediction of cluster constraints with a possible future 
mass bias prior of 1 - b = 0.80 + 0.01. The black shaded re¬ 
gion and dashed contours reproduce the curi'ent primary CMB 
and primary CMB-i-BAO constraints from Planck for the base 
ACDM model. The red shaded contours present the constraints 
expected from this mass bias prior applied to the present Planck 
cluster cosmology sample with the SZ-hBAO+BBN data set. 


meaning that the primary CMB has sufficient statistical weight 
to strongly ovenide the prior. The lensing power spectrum, in 
favouring slightly lower erg, reinforces the cluster trend so that a 
peak appears in the posterior for 2 niy in the red curve; it is not 
enough, however, to bring the posterior on the mass bias in line 
with the prior. This indicates that the tension between the cluster 
and primary CMB constraints is not fully resolved. 

One may then ask, how tight must the prior on the mass bias 
be to make a difference? To address this question, we performed 
an analysis assuming a projected tighter prior constraint on the 
mass bias. The informal target precision for cluster mass cali¬ 
bration with future large lensing surveys, such as Euclid and the 
Large Synoptic Survey Telescope, is 1%, and we consider the 
impact of a prior of 1 - I? = 0.80 + 0.01 on the present Planck 
cluster cosmology sample in Figs. 12 and 13. 

The latter figure compares the constraints from cluster counts 
for this mass bias to the present primary CMB constraints in the 
(Qm, crg)-plane for the base ACDM model. The bold, black dot¬ 
ted curve in Fig. 12 shows the predicted posterior on the neutrino 
mass from a joint analysis of the present Planck cluster counts 
and primary CMB. The same prior on a much larger catalogue 
would demonstrate a corresponding increase in sensitivity to 
neutrino mass. This simple projection highlights the importance 
and value of the more precise cluster cosmology expected in the 
future, and it provides clear motivation for significant effort in 
mass calibration. This effort will continue with larger samples 
of clusters with gravitational shear measurements, and also with 
the new technique of CMB lensing cluster mass measurements. 

7. Summary and discussion 

Our 2015 analysis incorporates a number of improvements and 
new information relative to our first study in Planck Collabora¬ 
tion XX (2014). With more data, we have a larger cluster cos¬ 
mology sample, increased by more than a factor of two, and 
we implement a two-dimensional likelihood over the counts in 


both redshift and signal-to-noise. We have also performed new 
tests of the selection function using MCXC and SPT cluster cat¬ 
alogs as truth tables. The selection function from these exter¬ 
nal checks and internal simulations of the Planck catalogue con¬ 
struction agree with each other and can be reasonably modelled 
by a simple analytical expression derived by assuming noise is 
the dominant factor (see the Appendix). One possible system¬ 
atic effect that warrants further study is IR emission from cluster 
member galaxies. Finally, we have examined the implications of 
three recent determinations of the cluster mass bias parameter, 
\ - b. The two-dimensional likelihood with the 2015 catalogue 
and mass bias priors will be implemented in CosmoMC. 

Our analysis confirms the results of the 2013 study. The 
counts are consistent with those of 2013, illustrated by the agree¬ 
ment in the constraints on Om and erg when using the same SZ 
observable-mass relations (see Fig. 3). The gain in statistical pre¬ 
cision is less than expected from the larger catalogue, which 
is likely related to the fact that the fit to the redshift distribu¬ 
tion with the X-ray prior on a is only marginal. Our new two- 
dimensional approach yields consistent, but more robust con¬ 
straints than the one-dimensional likelihood over just the red¬ 
shift distribution; it is less sensitive to the slope of the scaling 
relation, a, and it provides a better fit to the counts than in the 
one-dimensional case. 

Using the two-dimensional likelihood as our baseline, we ex¬ 
tracted new cosmological constraints using three different clus¬ 
ter mass scales represented by the mass bias prior distributions 
given in Table 2. The first two come from galaxy shear obser¬ 
vations of samples of Planck clusters. They differ by about Icr, 
with the WtG result favoring larger mass bias. We have also im¬ 
plemented a novel method for measuring cluster masses based 
lensing of the CMB temperature anisotropies behind clusters 
(Melin & Bartlett 2014). It gives a mass bias averaged over the 
entire cluster cosmology sample, although with larger statistical 
uncertainty. 

As a new method requiring further exploration, we con¬ 
sider CMB lensing less robust at present than galaxy lensing 
mass measurements, but highly promising. Similar CMB-based 
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mass measurements have recently been published by SPT (Bax¬ 
ter et al. 2014) and ACT (Madhavacheril et al. 2014). The ap¬ 
proach is appealing because it is subject to different systematic 
effects than gravitational shear and because it can be applied to 
large cluster samples thanks to the extensive sky coverage of the 
CMB experiments, with Planck of course covering the entire sky. 
Gravitational shear surveys will soon attain large sky coverage in 
the near future with the Dark Energy Survey (DES), and in the 
more distant future with the Euclid and WEIRST space missions 
and the Large Synoptic Survey Telescope. 

Our central result from analysis of the 2015 Planck cluster 
counts is shown in Eig. 7. Depending on the mass bias prior, 
we find varying degrees of tension with the primary CMB, as 
in 2013. The mass bias required to bring the cluster counts and 
CMB into full agreement is larger than indicated by any of the 
three priors and corresponds to 1 - = 0.58 + 0.04. Eig. 8 illus¬ 

trates the situation. The WtG prior almost eliminates the tension, 
but not quite, while both the CCCP and CMB lensing priors re¬ 
main in noticeable tension. Our largest source of modelling un¬ 
certain is, as in 2013, the mass bias. 

Tension between low redshift determinations of erg and the 
Planck primary CMB are not unique to the Planck cluster counts. 
Among SZ cluster surveys, both SPT and ACT are in broad 
agreement with our findings, the latter depending on which SZ- 
mass scaling relation is used, as detailed in our 2013 analysis 
Planck Collaboration XX (2014). And the new SPT cosmolog¬ 
ical analysis (Boequet et al. 2014) shows a significant shift be¬ 
tween the cluster mass scale determined from the velocity dis¬ 
persion or ix and what is needed to satisfy Planck or WMAP9 
CMB constraints (e.g.. Pig. 2 Boequet et al. 2014). In a study 
of the REPLEX X-ray luminosity function, Bohringer et al. 
(2014) also report general agreement with our cluster findings. 
On the other hand, Mantz et al. (2014) find the their X-ray cluster 
counts, when using the WtG mass calibration, match the primary 
CMB constraints. 

The situation is thus not yet satisfactory. It is unclear if these 
modest tensions arise from low-level systematics in the astro- 
physical studies, or are the first glimpse of something more im¬ 
portant. The most obvious extension to the base ACDM model 
that could in principle reconcile the differences is a non-minimal 
sum of neutrino masses. This, unfortunately, does not provide 
such a straightforward solution. While it is true that adding neu¬ 
trino mass does lower erg relative to the base ACDM prediction 
from the primary CMB, it does so at the the cost of increasing 
tension in other parameters; for example, it lowers Planck’s al¬ 
ready rather low value for the Hubble parameter. 

Pigure 14 highlights these points by showing constraints in 
the (Qm, CTg)- and the (Hq, crg)-planes for each of the mass bias 
parameters. Adding variable neutrino mass relaxes constraints 
from the primary CMB (shaded contours, which are the same in 
all three pairs of panels) towards lower erg, but by simultaneously 
increasing Qm and decreasing Hq. The remaining tension is mild 
in the case of the WtG mass prior, but more pronounced for the 
other two mass priors regardless of the neutrino mass. 

Another possibility is that baryonic physics influences the 
late-time evolution of the density perturbations. Strong feedback 
from active galactic nuclei (AGN) (Nagai et al. 2007; van Daalen 
et al. 2011; Martizzi et al. 2014) can potentially damp growth 
and lower erg through expulsion of matter from dark matter ha¬ 
los. This same effect could also reduce the mass of cluster halos 
and hence the prediction for their abundance, which is based on 
dark matter only simulations. It does not appear, however, that 
these effects are sufficiently large to explain the tension between 
low redshift and primary CMB constraints hinted at by the dif- 



Pig. 14: Effects of neutrino mass. The open magenta contours 
(solid, dashed and dotted) reproduce our cluster constraints with 
different data combinations for the base ACDM (i.e., Y^niy - 
0.06 eV). The violet shaded contours trace the constraints on the 
base ACDM model, while the other shaded regions give con¬ 
straints from the primary CMB combined with lensing and BAG 
when adding and marginalizing over variable neutrino mass. 


ferent observations. In addition, the violent feedback necessary 
for important impact might be difficult to reconcile with obser¬ 
vations of the baryon content of dark matter halos (e.g., Planck 
Collaboration Int. XI 2013). 

As conclusion, we return to the main uncertainty in interpre¬ 
tation of the cluster counts, namely the mass bias. It could be 
argued that the current accuracy is at the level of ~ 10 - 15%, 
based on the difference between different analyses and some¬ 
what larger than their quoted statistical uncertainties. Progress 
will certainly follow with improvement in these measurements. 
We illustrate the potential impact of a 1% determination of the 
mass bias in Pigs. 12 and 13. Such a result would, depending on 
the central value, significantly clarify the extent of any tensions 
and possible necessity for extensions to the base ACDM model. 
This precision is the avowed target of the large lensing surveys, 
such as Euclid, WEIRST and LSST. In the shorter term, we may 
expect valuable movement in this direction from DES and CMB 
lensing cluster mass measurements. 
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Fig. A.l: Robustness to the choice of cluster sample and the se¬ 
lection function model. The shaded contours give the cosmolog¬ 
ical constraints from the 2015 MMF3 cluster catalogue using 
the analytical (Error Function, ERF) selection function model 
(grey), the MMF3 Monte Carlo selection function (red), and 
the Monte Carlo selection function for the intersection sample 
(blue). Our final constraints are obtained from the MMF3 ERF 
model. For this comparison, we adopt the SZ-i-BAO-i-BBN data 
set and we fix (1 - /?) = 0.8. 



Appendix A: Modelling uncertainties 

We examine the robustness of our cosmological constraints to 
modelling uncertainties. We first consider sensitivity to the cos¬ 
mological sample and to our modelling of the completeness 
function in Sect. A.l, and then look at the effect of using a dif¬ 
ferent mass function in Sec. A.2. In Sec. A.3, we show that our 
constraints are robust against redshift evolution of the scaling 
relations. 

Appendix A. 1: Choice of the sampie and seiection function 

For our baseline analysis, we use the MMF3 cosmological sam¬ 
ple and its associated completeness based on the analytical ap¬ 
proximation using the Error Function (Eq. 14). In Fig. A.l 
we show how the Monte Carlo determined selection function 
changes the cosmological constraints (labelled QA for "Quality 
Assessment" in the figure). We also present the constraints ob¬ 
tained from the intersection sample defined in Sec. 2. The figure 
is based on the ID N(z) likelihood, for which the Monte Carlo 
completeness can be easily computed, and we use the baseline 
SZ-hBAOh-BBN data set and fix (I - b) ^ 0.8. The MMF3 ERF 
contour is thus close to the q - 6 contour of Fig. 3, the only 
difference being that crin y is fixed to zero in Fig. A. 1. while it is 
constrained by the Table 1 prior in Fig. 3. The impact of adopting 
the intersection sample and/or the QA completeness function is 
small (< 0.5 (t) for both Qm and erg. 

Appendix A.2: Mass function 

We use the Tinker et al. mass function for our baseline anal¬ 
ysis. To characterize the influence of this choice, we examine 
constraints when adopting the mass function from Watson et al. 
(2013) instead. We employ our 2D N{z,q) likelihood (with the 
CCCP mass bias prior and a constrained) and combine with 
BAO and BBN prior constraints, and show the result in Fig. A.2. 
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Fig. A.2: Robustness to the choice of mass function. The grey 
shaded contours give the cosmological constraints when using 
the Tinker et al. mass function, corresponding to our final result. 
This is compared to constraints obtained when using the Watson 
et al. mass function, shown as the red shaded contours. In this 
figure we adopt the SZ-hBAOh-BBN data set and the CCCP mass 
bias prior. 


The Tinker et al. contour of Fig. A.2 is thus identical to the 
N{z,q) contour with a free, as given in Fig. 5. The new mass 
function shifts our constraints by about Icr towards higher Qn, 
and lower erg, along the main degeneracy line, hence increas¬ 
ing the tension with the Planck primary CMB constraints. Note 
that we use the general fit from Eq. (12) of Watson et al. (2013) 
(independent of redshift). This was not the case for our 2013 pa¬ 
per (Planck Collaboration XX 2014) where we adopted the AHF 
fit with parameters varying with redshift in the first ArXiv ver¬ 
sion of the paper, which was subsequently found to be incorrect 
(fourth version of the paper on ArXiv). 
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Fig. A.3: Robustness to redshift evolution in the SZ-tnass scal¬ 
ing relation. The different contours show the constraints when 
relaxing the redshift evolution exponent,/?, of Eq. (7). The black 
contours result from hxing jS - 0.66, our hducial value through¬ 
out, with a constrained by the Gaussian X-ray prior of Tab. 1 . 
Applying a Gaussian the prior on /? instead, from Tab. 1, pro¬ 
duces the blue contours, while the red contours result when we 
also leave a free. In this figure we adopt the SZh-BAO-i-BBN 
data set and the CCCP mass bias prior. 


Appendix A.3: Redshift evoiution of the Y-M reiation 

Throughout our baseline analysis, we fix the redshift evolution 
exponent fi - 0.66 (self-similar prediction) in Eq. (7). Here we 
examine the impact of allowing this parameter to vary. Con¬ 
straints when leaving p free are shown in Fig. A.3. The "P hxed, 
a constrained" case corresponds to the 2D N{z,q) likelihood 
(CCCP mass bias prior and a constrained) combined with BAG 
and BBN, as in Fig. A.2. This contour is also identical to the 
N{z,q) contour with a free, shown in Fig. 5. For the "P con¬ 
strained" cases, p is allowed to vary over the range 0.66 ± 0.50 
(Table 1). This increases the size of our constraints along the ma¬ 
jor degeneracy between Q,„ and erg, but does not bring them into 
any closer agreement with the primary CMB. 
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